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Abstract 

The nuclear collective response at finite temperature is investigated 
for the first time in the quantum framework of the small amplitude 
limit of the extended TDHF approach, including a non-Markovian col- 
lision term. It is shown that the collision width satisfies a secular equa- 
tion. By employing a Skyrme force, the isoscalar monopole, isovector 
dipole and isoscalar quadrupole excitations in °Ca are calculated and 
important quantum features are pointed out. The collisional damping 
due to decay into incoherent 2p-2h states is small at low temperatures 
but increases rapidly at higher temperatures. 
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1 Introduction 



After discovery of giant dipole resonance in 1947, much work has been done 
to understand the properties of nuclear collective vibrations built on the 
ground state and excited states. Most of these theoretical investigations are 
based on the RPA theory which is quite successful in describing the mean 
resonance energies and fragmentation of the excitation strengths at zero and 
finite temperatures. However, the RPA approach, which is in fact the small 
amplitude limit of the TDHF theory, is not suitable for describing damping 
of the collective excitations |1J. Damping arises mostly by mixing of the 
collective state with the nearby complex states ||. The importance of the 
intrinsic compound nucleus life-time has also been stressed ||, |J]. As a result 
of the mixing with complex states, the excitation strength spread around 
the mean resonance energy, and furthermore the damping width increases 
with the intrinsic temperature of the system as observed in giant dipole res- 
onances in 120 Sn [[| and 208 Pb nuclei |6j, [7|, ||. In order to describe the 
nuclear collective response including damping, it is necessary to go beyond 
the RPA theory by incorporating coupling between the collective states and 
the doorway configurations. There are essentially two different approaches 
for this purpose: (i) Coherent mechanism due to coupling with low-lying sur- 
face modes which provides an important mechanism for damping of giant 
resonances in particular at low temperatures || £|, (ii) Damping due to the 
coupling with incoherent 2p-2h states which is usually referred to as the col- 



lisional damping ||10|| . The small amplitude limit of the extended TDHF is 
an appropriate basis for investigating collective response, in which damping 
due the incoherent 2p-2h decay is included in the form of a non-Markovian 
collision term |TT], [T2], [13]]. Based on this approach, the incoherent contribu- 
tion to damping at finite temperature has been calculated in Thomas-Fermi 
approximation in refs. |TJ, |Tj], IB . Calculations using the Markovian limit 



of this semi-classical treatment, the so-called Boltzmann-Uehling-Uhlenbeck 
(BUU) approach, are discussed by many authors (for a review see JHJ). How- 
ever, as far as the collective behavior of nuclei at moderate temperature are 
concerned one may worry about the adequacy of semi-classical calculations 
which neglect most of the quantum features but the Pauli principle. 

In this work, we present a first quantal investigation of the nuclear collec- 
tive response at zero and finite temperatures on the extended TDHF frame- 
work in small amplitude limit, which may be referred as an extended RPA 
approach. In this approach, in contrast to the semi-classical treatments, the 



2 



shell effects are incorporated into the strength distributions as well as the col- 
lisional damping widths. We point out that the damping widths should be 
calculated by solving a secular equation. We compute the isoscalar monopole, 
isocalar quadrupole and isovector dipole strength distributions in m Ca at fi- 
nite temperatures by employing an effective Skyrme force. 



2 Collective response at finite temperature 

In the extended TDHF theory, the evolution of the single particle density 



matrix pit) is determined by a transport equation [16, 18, 19, |2C 

El, 

ih^p - [h{p),p] = j* drTr 2 [v, G(t, t - r)F 12 (t - r)G\t, t - r)] (1) 

where h(p) is the mean-field Hamiltonian and the right hand side represents 
a non-Markovian collision term with 
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(1 - pi)(l - P2)vpip 2 - Plp2v(l ~ Pl)(l - p 2 ) (2) 



and G(t, t — r) = T ■ exp[— 1 f t _ r dt'h(t')} denotes the mean-field propagator. 
The small amplitude limit of the extended TDHF theory provides a suit- 
able framework to describe collective vibrations including damping due to 
the coupling with incoherent 2p-2h excitations. The small deviations of the 
single-particle density matrix 5 pit) = pit) — po around a finite temperature 
equilibrium state p is determined by 

iU^Sp _ [ ho: S p] - [5U + F, po] = I • Sp (3) 

where 5U = (dU/dp)o ■ 5p represents small deviations in the effective mean- 
field potential and F(r,t) = F(r) exp(—iut) + h.c. is a one-body excitation 
operator with harmonic time dependence. An explicit expression of the lin- 
earized form I ■ 5p of the non-Markovian collision term can be found in a 
recent publication [|TB| . 



The linear response of the system to the external perturbation F is deter- 
mined by expanding the small deviation Sp(t) in terms of finite temperature 
RPA modes, 

Sp(t) = $> A (t)pl + zl(t) Px } (4) 

A>0 
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where the finite temperature RPA modes, p\ and p x , are defined by 

hu x p\ - [ho, p[] - [h[, po] = 0. (5) 

Here to x is the mean-frequency of the RPA mode and h\ = (dU/dp) ■ p\ 
represents the positive frequency part of the vibrating mean-field. It is con- 
venient to introduce the collective operators, O x and O x , associated with the 
RPA modes as p\ = [0\, p ] and p x = — [O x ,p Q ], and they are orthonormal- 
ized according to Tr[O x ,Oj l ]p = Substituting the expansion (4) and 
projecting by O x , we find that the amplitudes of the RPA modes execute 
forced harmonic motion, 

- ihjzx + (Hujx - t -T x )z x =< [O x , F] > (6) 

where < [0\,F] > = Tr[0\, F]p and T x denotes the collisional damping 
width of the mode. Here, we neglect a small shift of the mean frequency 
uo\ arising from the principle value part of the collision term. Solving this 
equation by Fourier transform, the response of the system to the external 
perturbation F can be expressed as 

8p(uj) — R(uj, T) ■ F (7) 

where R(u,T) denotes the finite temperature extended RPA response func- 
tion including damping 

d i, , T \ - i < i\P\\j >< k \px\ l > < k \p\\ l >< i\P\\j > \ /on 
R l3M {u, T) - ^ \- huj _ hux + Fx + nu + hux + , Vx ) ■ (8) 

The strength distribution of the RPA response is obtained by the imaginary 
part of the response function, 

S(u,T) = --Tr{F^ImR{uj,T)F} (9) 

7T 

= "E{l < >o \ 2 D(oo-u; x )-\< [0{,F] > \ 2 D{uj + u x )} 

71 \>o 

where the sum goes over the positive frequency modes and 

^-^' CTiiW ' (10) 
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The main features of the strength function are usually discussed in terms of 
sum rules, which are calculated from the RPA response as, 



oo 

k 



m k (T) = / uj h dujS(uJ,T) (11) 
Jo 

However, due to the Lorentzian shape of the poles these moments are not well 
defined for k > 1. For a Hermitian excitation operator the energy weighted 
sum rule for k — 1 is not effected by the damping and it is given by, 

mi (T) = J2^\<[Ox,F]> \ 2 . (12) 

A>0 

For a multipole operator, the energy weighted sum rule leads to m\(T) = 
| < [F*, [H,F]] >o which is exactly satisfied by the finite temperature RPA 



sum rule, as shown by Vautherin and Vinh Mau 25 



In the Hartree-Fock basis the finite temperature RPA equation reads, 

(huox-ti + tj) <i\0{\j>+ (13) 
< ik\v\jl > A (n, - n k ) < l\0{\k >= 

l^k 

where v = (dU/dp)o, the indices i,j,.. represent all single particle quantum 
numbers including spin-isospin, and n k = 1/[1 + exp(efc — ep)/T\ denotes 
the finite temperature Fermi-Dirac occupation numbers of the Hartree-Fock 
states. At zero temperature these occupation numbers are zero or one, so 
that the RPA operators O x , 0\ have only particle-hole and hole-particle 
matrix elements. At finite temperatures the RPA functions involve more 
configurations including particle-particle and hole-hole states. By associating 
a single index with the pair of indices (i,j), the RPA functions can be regarded 
as a vector, and in this manner eq.(13) can be expressed as an eigenvalue 
equation for finite temperature RPA modes p3| . According to the small 



amplitude limit of the extended TDHF equation, the damping width of RPA 
modes due to decay into incoherent 2p-2h doorway excitations is given by 



16 



1 H I < ij\[0\,v\\kl > A \ 2 Di jM [n h ninirij - nin^ni] (14) 
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where Hi = 1 — n,. In reference |]16| , neglecting the damping of the collective 
amplitude in the collision term, the energy conserving factor is taken as a 
sharp delta function as D^^x = Im{huj\ — — €j + e k + e\ — ir})~ x Here, we 
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take into account for depletion of the collective amplitude in the collision 
term by substituting u\ — |T A in place of oj\. Then, the factor takes a more 
appropriate Lorentzian form, 

DiiM = (^ A -e l - ej +e fc / +Q) 2 + (r A /2)2 (15) 

Then, expression (14) becomes a secular equation for the damping width. As 
we will see this self-consistency is of major importance in order to properly 
compute the collision width. The collective mode damps out by mixing with 
the intrinsic states of increasing complexity. The sequence of the complexity 
of the states can be classified according to the exciton number as mixing 
with 2p-2h, 3p-3h,- • •, Np-Nh,- • • states. The expression (14) contains only 
the mixing with 2p-2h doorway states in accordance with the extended TDHF 
theory. In order to incorporate the damping due to mixing with more com- 
plex states, the extended TDHF should be improve by including higher order 
correlations beyond binary collision term. The effect of the higher order mix- 
ing may be approximately taken into account by introducing an appropriate 
decay width, of 2p-2h states in the expression (15) 

ij ' kl = (hu x - Ci - 7 j + e k + e ,) 2 + «r A + r y , fcI )/2) 2 (16) 

Then, the secular equation can also be solved considering these higher order 
effects. 



3 Results 

We calculate the isoscalar monopole, isoscalar quadrupole and isovector dipole 
excitations in m Ca at several temperatures. We use the Skyrme interac- 



tion SGII for the Hartree-Fock and RPA calculations |26] and we neglect the 
temperature dependence of single particle energies and wave functions. We 
determine the hole states by solving the Hartree-Fock problem in coordinate 
representation. Then, the particle states are generated by diagonalizing the 
Hartree-Fock Hamiltonian in a large harmonic oscillator representation by 
including 12 major shells. In this manner, unbound continuum states are 
approximately included in the RPA calculations. The RPA strength distri- 
butions of the monopole F (r) = r 2 , dipole i*i(r) = r z rYi (r) (in isospin 
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symmetric systems N = Z), and quadrupole F 2 (r) = r 2 Y" 2 o(r) excitation 
operators at temperatures T = 0,2,4 MeV are shown in figure 1. As seen 
from the top panel of figure 1, the monopole strength at T = MeV ex- 
hibit a large Landau spreading over a broad energy region E = 16 — 28 MeV 
with an average energy E = 21.5 MeV. The recent experimental data also 
show a broad resonance around a peak value of 17.5 MeV|27||. For increasing 
temperature, transition strength spread a broader range towards lower en- 
ergies. As shown in the middle panel, the strength distribution of isovector 
dipole shows a weaker temperature dependence than monopole. At T = 0, 
the dipole strength is concentrated at range E = 16 — 23 MeV. The Landau 
width is large and is spreading for increasing temperature. However, the 
average energy of the main peak remains nearly constant around E = 16.5 
MeV. The experimental data shows a broad resonance at around 20 MeV 
with a width close to 6 MeV. As illustrated at the bottom panel of figure 



1, the RPA result at T = MeV gives a very collective quadrupole mode 
peaked at E = 17.5 MeV, which agrees well the experimental finding of an 



average energy 17 MeV [29| and the calculations of Sagawa and Bertsch [30 



At higher temperatures in addition to p-h excitations, p-p and h-h excita- 
tions become possible. The p-p and h-h configurations mainly change the 
strength distribution at low energy side at E = 4 MeV. As a result, the giant 
resonance has less transition strength. 

Figure 2 illustrates the average energy < E >= mi/m (solid lines), and 
the ratios (mi/m.i) 1 ' 2 (short dashed lines) and {m^,/mi) 1 ^ 2 (long dashed 
lines) of the monopole, dipole and quadrupole excitations as a function of 
temperature. The behavior of the average energy of the monopole resonance 
is particularly interesting, since it may be related to the compressibility co- 
efficient of nuclear matter [pi], |32]j . 

We obtain the collisional damping widths of the collective states by cal- 
culating the expression (14) and solving the associated secular equation by 
graphical method. We note that the sums over single particle states needed 
to evaluate the expresion (14), have been performed explicitely using the pro- 
jection of the total spin, m, as one of the explicit quantum numbers as done 
in ref. ||33|| . Figure 3 illustrates examples of the graphical solution for giant 
dipole and quadrupole excitations at two temperatures T = 0, 2 MeV. In this 
figure, the curves with dashed lines are obtained by calculating the right hand 
side of (14) as a function of T x = T in . The intersection of this curve with 
the diagonal line determines the solution. The effect of the damping width 
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^ij,ki of 2p-2h states may be approximately incorporated by taking a larger 
value of T in as indicated in eq.(16). As seen from figure 3, in most cases, the 
selfconsistent value of the damping width saturates very rapidly, and hence 
it is not modified very much by increasing Ti n . Figure 4 shows the damping 
widths as a function of temperature, that are averaged over several nearby 
states with strengths more than 10% of the EWSR. The results for monopole, 
dipole and quadrupole are indicated in the top, middle and bottom panels, 
respectively. Collisional damping widths are generally small at low temper- 
atures, but rapidly grow for increasing temperature. This increase appears 
to be more complex than the semi-classical quadratic prediction. Depending 
upon the mode, the increase may be linear or may saturate. 

In order to understand this behavior, it is convenient to write the ex- 
pression (14) of the damping width as sum over energy bins in energy E = 
e i + e j ~~ e k ~ e i °f 2p-2h states, 

rx = W92 P -2 h (E)W x (E)D(hcu x -E) (17) 

1 E 

where each bin has a small energy interval AE around E, and D(huo\ — E) 
is the Lorentzian factor given by eq.(15). Here W\(E) denotes the average 
transition rate, 

W *( E ) = -TFtXH < U'|[°A,upJ >a f\n k n l fi i fi j - mnjUkfii]. (18) 

g2p-2h{ti) AE 

In this expression many terms are vanishing either due to the selection rules or 
due to the Pauli blocking factors. The quantity g 2p -2h(E) is the total number 
of 2p-2h states in the energy interval including only those states which are not 
Pauli blocked and which have non- vanishing matrix elements of the transition 
rate, i.e., which can be coupled to the phonon quantum numbers. Figures 
5 and 6 show the density of 2p-2h states and the average transition rates 
as a function of the 2p-2h energy for dipole and quadrupole excitations at 
T = MeV and T = 3 MeV. In dipole mode, there is no odd parity 2p-2h 
states available in the vicinity of the collective energy, hence the average 
transition rate W\{E) and the density of states vanish at zero temperature. 
As a result, the collisional damping of giant dipole in m Ca is zero at T = 
MeV. This behavior is a particular quantum feature due to shell effects in 
the extended RPA calculation of double magic light nuclei, and it can not be 
described in the framework of semi-classical approaches. In medium weight 
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and heavy nuclei, in the vicinity of the GDR strength, there are few odd 
parity 2p-2h configurations involving intruder states associated with the spin- 
orbit coupling. As a result, we expect to find a small finite damping of 
giant dipole resonance at zero temperature. For increasing temperature, the 
available phase space becomes much larger and the collisional damping of 
the GMR and GDR increases. This is not the case for the GQR because the 
increase of the phase space is compensated by a reduction of the magnitude of 
average transition rates. As a result, the damping width of giant quadrupole 
appears to saturate above T = 3 — 4 MeV. 

Figure 7 shows the strength distributions including the collisional damp- 
ing. The giant dipole strength at T = MeV is smoothed by performing an 
averaging with a Lorentzian weight with a width of 0.5 MeV. The excitation 
strengths become broader for increasing temperature. The peak position of 
the monopole resonance does not change much, but the peak position of 
dipole slightly shifts down and quadrupole slightly shifts up in energy. This 
is a signature of the reduction of the collectivity of those states with tem- 
perature because the peak energy moves back towards the single particle 
expectations. 

Figure 8 illustrates the FWHMs of the collective excitations at several 
temperatures. The widths of giant monopole and giant dipole at T = MeV 
are mainly due to fragmentation of the collective strength, i.e. the so-called 
Landau spreading, which is about 4 MeV in both cases. Since it is difficult 
to extract well defined values, the FWHM of these modes at T = MeV are 
left open in the figure. The total widths increase further by mixing of the 
collective mode with incoherent 2p-2h excitations at higher temperatures. 
However, the total width does not present the parabolic behavior predicted 
by semi-classical calculations. 

In figure 9, long-dashed lines and solid lines show the integrated strengths 
over the energy interval 10 — 40 MeV as a function of temperature in the 
RPA and the extended RPA, respectively. As a reference, the total strength 
mi is also indicated by short-dashed lines. In the RPA calculations, the 
modes retain a high degree of collectivity even at temperatures at T = 4 — 
5 MeV. However, in the extended RPA approach, as a result of damping, 
the excitation strength become broader and the collectivity diminishes for 
increasing temperature. 
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4 Conclusions 



We investigate isoscalar monopole, isoscalar quadrupole and isovector dipole 
excitations of 40 Ca at finite temperature in the basis of the small amplitude 
limit of the extended TDHF. The extended TDHF goes beyond the thermal 
RPA approach by including damping due to decay into incoherent 2p-2h exci- 
tations. We calculate the excitation strength distributions in a self-consistent 
Hartree-Fock representation by employing a Skyrme force with SGII param- 
eters. At T = 0, the monopole and dipole strengths are fragmented and 
spread over a broad range because of the Landau damping, while quadrupole 
exhibits a single peak structure. For increasing temperature, strength in all 
cases becomes broader and hence the collectivity is reduced. The incoherent 
damping widths at low temperatures are, in general, small and thus leaving 
room for a possible coherence effect of doorway states in the description of 
the damping properties. At high temperature the collisional damping be- 
comes large and may even dominate the spreading width since the coherence 
effect is expected to diminish rather rapidly. For increasing temperature, the 
collisional damping predicted by the quantal calculations, evolves in a more 
complex manner than the quadratic increase predicted by the semi-classical 
calculations. An interesting property of the collisional damping is that it 
may saturate for increasing temperature. In fact, our calculations indicate 
that the damping width of giant quadrupole saturates around T = 3 — 4 
MeV, however a saturation of the giant monopole and dipole modes is not 
visible at these temperatures. There are important quantal effects in the 
collective behavior of a hot nuclear system as illustrated in P3 |. Investi- 
gations presented here, also indicates that, the quantal effects has a large 
influence on the damping properties of collective excitations at low temper- 
atures, which may even persist at relatively high excitations. As illustrated 
in ref. |T^], the magnitude of the collisional damping is rather sensitive to 
the effective residual interactions, for which an accurate information is not 
available. The effective Skyrme force is well fitted to describe the nuclear 
mean-field properties, but not the in-medium cross-sections and damping 
properties. Therefore, a systematic study of the effective interactions in this 
context is clearly called for. However, our investigation, while remains semi- 
quantitative, gives a valuable insight on the quantal properties of collective 
excitations at finite temperature. 
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Figure 1: RPA strength distributions in Ca 40 as a function of the energy 
at temperatures T = 0,2,4 MeV for isoscalar monopole + (top), isovector 
dipole 1~ (middle) and isoscalar quadrupole 2 + (bottom) excitations. 
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Figure 2: Averaged energy of the monopole, dipole and quadrupole excita- 
tions in 40 Ca mi/rriQ (solid line) and the moment ratios (mi/m_i) 1 '' 2 (short 
dashed line) and (m 3 /mi) 1 ^ 2 (long dashed line) as a function of the temper- 
ature. 




Figure 3: Graphical solution of the secular equation for the damping width 
for L n = 1" (top) L n = 2+ (bottom) at T=0 MeV (left) and T=2 MeV 
(right). 
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Figure 4: Collisional damping widths that are averaged over nearby states 
with more than 10% of the EWSR,lft>r monopole, dipole and quadrupole 
modes as a function of temperature. 




10 20 30 40 50 10 20 30 40 50 



E (MeV) E (MeV) 



Figure 5: Top: Energy dependence of the density of 2p-2h states g2 P -2h for 
L 7r =l~(left) and L 7T =2 + (right) at zero temperature. Bottom: Averaged 
transition rate between collective states and 2p-2h states as a function of the 
energy of 2p — 2h states. 
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Figure 6: The same as in figure 5, but at T = 3 MeV. 
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Figure 7: The extended RPA strength distributions of the monopole, dipole 
and quadrupole excitations as a function of temperature. 
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Figure 8: The full widths at half maximum of the strength distributions at 
several temperatures. 
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Figure 9: Integrated strengths over the energy interval 10 — 40 MeV in 
the RPA (long-dashed line) and the extended RPA (solid line). The total 
strengths mi are plotted as a referenSfi by short dashed lines. 



